Insulator-to-metal phase transition in Yb-based Kondo insulators 
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The periodic Anderson lattice model for the crystalline electric field (CEF) split 4f quartet states 
is used to describe the Yb-based Kondo insulators/semiconductors. In the slave-boson mean-field 
approximation, we derive the hybridized quasiparticle bands, and find that decreasing the hybridiza- 
tion difference of the two CEF quartets may induce an insulator-to-metal phase transition. The 
resulting metallic phase has a hole and an electron Fermi pockets. Such a phase transition may 
be realized experimentally by applying pressure, reducing the difference in hybridization of the two 
CEF quartets. 

PACS numbers: 71.27.+a, 75.30.Mb, 75.20.Hr 



Kondo insulators or semiconductors, such as YbBi2, 
belong to strongly correlated electron systems 1 -^, in 
which the conduction electrons hybridize with the local- 
ized 4f-electrons and the strong Coulomb repulsion re- 
sults in highly renormalized quasiparticle bands with a 
small indirect energy gap 3-6 . To study the characteristic 
properties of these materials at low temperatures, some 
experiments have been attempted to make the insulat- 
ing gap vanish by applying an external magnetic field 7 ' 8 
or pressure^, leading to an insulator-metal phase transi- 
tion. Such a transition under the external magnetic field 
has been considered in the previous studies 1 ^ - — , how- 
ever, the microscopic mechanism for the pressure induced 
insulator-to-metal transition remains far from being fully 
understood. 

For the Yb-based Kondo insulators, an external pres- 
sure can affect the hybridization between 5d band elec- 
trons and the more atomic-like 4f electrons, giving rise to 
the intermediate valence behavior. The Yb valence is di- 
rectly related to the number of 4f- holes by v = 2 + rih- 
At the ambient pressure, i%h spans a broad range between 
and 1 in Yb-based compounds, and the intermediate va- 
lence reflects the hybridization of the energetically close 
Yb 2+ (4f 14 ) and Yb 3+ (4f 13 ) configurations. The elec- 
tronic configuration of Yb 3+ (4f 13 ) can be regarded as a 
single hole in the 4f-shell, while the configuration of Yb 2+ 
(4f 14 ) corresponds to the closed 4f-shell. Taking into ac- 
count the much larger strength (1.3eV) of the spin-orbit 
coupling 14 , a j = 7/2 f-hole state is split into a quar- 
tet and two doublet states by the crystalline electric field 
(CEF) under the cubic symmetry, which is the usual lat- 
tice structure of YbBi 2 . These two doublets are almost 
degenerate and may be treated as a quasi-quartet. Thus, 
a periodic Anderson lattice model with U — > oo for the 
CEF split 4f states can be used to describe these Yb- 
based Kondo insulators or semiconductors^ 5 .. 

It has been further pointed out that the anisotropic hy- 
bridizations of the two CEF quartets play an important 
role in the formation of the two dispersive spin resonances 



at the continuum threshold^ 5 -, the most salient features 
observed by inelastic neutron scattering experiments in 
YbBi2 (Ref|5||6|)- Motivated by this analysis, we further 
notice that, above a threshold of the CEF splitting, de- 
creasing the difference in the hybridization of the two 
CEF quartets may cause an overlap between the middle 
lower and upper hybridized quasiparticle bands, leading 
to an insulator-to-metal phase transition. Experimen- 
tally, this phase transition can be realized by applying 
pressure, reducing the difference in the hybridization of 
the two CEF quasi-quartets. Such a pressure induced 
insulator-to-metal transition has been observed in the 
Kondo insulator SmBg, where the electrical resistivity 
has been measured below 80 K and under pressure be- 
tween 1 bar and 70 kbar (Ref£) . Above the critical pres- 
sure 40 kbar, a transition occurs from a Kondo insulator 
to a metallic heavy fermion liquid and a non- Fermi liquid 
behavior has been found^. 

In this paper, we will carefully study the periodic An- 
derson model with U — > oo for the CEF split 4f quartet 
states. Using the slave-boson mean-field approximation, 
we will derive four quasiparticle bands resulting from the 
hybridization between the conduction electrons and lo- 
calized 4f-hole states, and find that decreasing the hy- 
bridization difference of the two CEF quartets indeed 
induce an insulator-to-metal phase transition. The re- 
sulting metallic phase has a hole and an electron Fermi 
pockets. By including the Coulomb interaction between 
the localized and conduction electrons, we discuss the 
possible instability of the resulting metallic phase. 

To describe the Yb-based Kondo insulators or semi- 
conductors, the periodic Anderson lattice model with 
U — > oo for the CEF split 4f states has been introduced 15 
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where the first term denotes the conduction electron 
band, the second term stands for the binding energy 
of the 4f-hole, and A 7 (Ai = 0,A 2 = A) is the 
CEF splitting energy for the two quasi-quartets with 
7 = (r, m), where T = 1,2 denotes the quartets and 
m = 1 — 4 represents the four-fold orbital degeneracy. 
Due to the exclusion of the double occupancy, a pro- 
jection has been implemented by using the slave-boson 
representation^. Then the Yb 2+ (4f 14 ) configuration 
without a 4f-hole state can be accounted for by an aux- 
iliary boson state b\ |0), while the Yb 3+ (4f 13 ) configu- 
ration with a 4f-hole state is represented by a fermion 
state fj |0). The conduction electrons hybridize with 
the /-hole at each lattice site in both quartets with 
different strengths. At each lattice site the constraint 
Qi = b\bi + J^-y fj-yfif = 1 nas to be enforced, and the 
total Hamiltonian is T~L + J2i ^i(Qi ~ 1): where Ai is the 
Lagrange multiplier. 

Now the slave-boson mean-field approximation is per- 
formed by neglecting the fluctuation of the Bose field 
(b\) = (bi) = b and the site dependence of the local field 
A, = A. Within these approximations, the mean-field 
model Hamiltonian can be written as 

H m f = ^[ekdk 7 rfk 7 +£ 7 /k 7 /k7 + ^7( d k 7 /k 7 + /i-C.)]+A/'eo, 

where e 7 = ef + A 7 + A is the renormalized energy level 
of the localized states, V 1 — bV 1 , and eo = A(6 2 — f). It 
should be noticed that the dependence of the hybridiza- 
tion strength on k has been neglected, i.e., Vk 7 = V-y 
Furthermore, we will replace V 1 by Vr for simplicity. By 
performing the Bogoliubov transformation 

C% 7 = ^k-yC^ + ^k/k 7 , Ac 7 = — ^k 7 ^k 7 + /^k 7 /k 7 , (3) 

we can diagonalize the quadratic Hamiltonian and obtain 
-Hmf = £ (E+ 7 at 7 a kl + E^^) , (4) 

k,7 

with four hybridized quasiparticle bands are 



where O(E^) is the step function. The chemical potential 
/i and the Lagrange multiplier A have to be determined 
self-consistently according to the conservation of the total 
number of particle per lattice site n c + rif — 2. Depend- 
ing on the parameter values £/, A, andVr, the variational 
parameters b and A are also determined self-consistcntly. 
From the hybridized quasiparticle band structure, the 
ground state of the system can be an insulating state, 
where the two lower bands are filled completely, leav- 
ing an indirect energy gap. As the k dependence in E g 
appears through the conduction electron energy e^, sum- 
mations over k can be transformed into an integral over 
energy e in the interval [— D,D]. By assuming a con- 
stant density of states, the ground-state energy is thus 
evaluated as 
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where Ap(x) 



(x ± £t) 2 + 4Vp 2 , e 2 = £i + A, and 

V2 = bV2 = b(V\ + 5V). Minimizing the ground-state 
energy density with respect to b and A, respectively, we 
obtain the following self-consistent equations 
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However, we notice that there exists another possible 
structure of the quasiparticle bands, where the chemi- 
cal potential 11 cuts through the two middle hybridized 
quasiparticle bands E^ x and E^ 2 at £1 and £21 respec- 
tively. Both these energy parameters are determined by 
the equation E^ = E^ 2 — fx. From the condition of the 
total number of particles per lattice site n c + rif — 2, we 
can derive the result £1 = —£2 = —£, and 
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while the Bogoliubov parameters /Xk 7 and i^k 7 are given 
by 



^k 7 
^k 7 



1 
71 



1 ± 



e k 



( e k - £ 7 ) 2 



4V 7 2 



1/2 



(6) 



These two parameters describe the contributions of the 
conduction electron band and localized /-hole band to 
the hybridized quasiparticles, respectively. 

Moreover, the ground-state energy per site is given by 



£ S = ^ [^(£k 7 ) + e^9{eq 
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2£ + A = A+(0+A 2 -(£). 



(10) 



Here £ can be used to characterize the insulator-to-metal 
transition. When < £ < D, the ground state should be 
metallic, while for £, = D the ground state corresponds to 
a critical point. The corresponding ground-state energy 
density in the metallic phase is thus expressed as 

E g = ^L[(3D - Oei + (D + Oea + £ 2 - D 2 ] 



V2 ln .A2 (0+^-^2 



2D A^(D) + D-e 1 2D A+(D) - D - e 2 
-L[(e + f 1 )A+(0 + (£'-£i)Ar(^)] 
g^[(C ~ e a )Aj(0 + (D + e 2 )A+(D)} + e . (11) 
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FIG. 1: (color online). The quasiparticle band structure for 
Vi = 0.4D, e/ = -0.573, and A = 0.173. (a) Insulator phase, 
(b) Critical point, (c) Metallic phase. 



By minimizing E g with respect to b and A, the corre- 
sponding self-consistent equations can be deduced to 
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(12) 

In order to deduce the ground state phase diagram, 
we should first numerically solve Eq.® for the insu- 
lating phase and Eqs. flTUl) and (IT21 for the metallic 
phase, respectively. The hybridized quasiparticle band 
energy versus the momentum along the diagonal direc- 
tion T (0,0,0)- > M (7r,7r,7r) are plotted in Figfflwith 
Vi = 0.4D, e/ = -0.573, and A = 0.173 for three differ- 
ent values of SV. As shown in FigQJa) f° r $V — 0.1873, 
there opens an indirect gap between the middle upper 
and lower bands, corresponding to an insulating phase. 
In FigQJb) for SV — 0.12673, the middle upper and lower 
bands just meet at the chemical potential, correspond- 
ing the critical point of the transition. Since we have 
£ = 73 at the critical point, the ground-state energies 
of the metallic and insulating phases are equal. So the 
insulator-metal transition is a continuous second-order 
phase transition. Finally, in FigOJc) for SV = 0.0173, 
the middle lower and upper bands overlap, and the chem- 
ical potential cuts through these two bands, which cor- 
responds to the metallic phase. 

The critical condition under which the insulator-metal 
transition occurs can be determined from Eq. (|10p and 
Eq.(fl~2"j) by setting £ = 73. Then the ground-state phase 
diagram can constructed for V\ = 0.473 and e/ = — 0.573 
and is shown in Figl^a). Clearly there exists a threshold 
of the CEF splitting energy A c , and only when A > A c 
the insulator-to-metal phase transition occurs by turn- 
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FIG. 2: (a) Ground state phase diagram for V\ = 0.473 and 
e/ = —0.573; (b) Indirect energy gap as a function of the 
hybridization difference of the two CEF quartets SV for A = 
0.173. 



ing the difference in hybridization of the two CEF quasi- 
quartets. The change of the indirect gap is another evi- 
dence to characterize the insulator-to-metal phase transi- 
tion, and can be also calculated and displayed in FigJ^b) 
for A = 0.173. It shows that the indirect quasiparti- 
cle gap decreases almost linearly with decreasing the hy- 
bridization difference of the two CEF quartets, and this 
energy gap finally vanishes at 5V C . There is another criti- 
cal value SV* , where the top energy levels of the two lower 
quasiparticle bands interchange with each other around 
the Brillouin zone boundary. Then the indirect energy 
gap has a cusp. 

Actually, such an insulator-to-metal phase transition 
can be realized experimentally. There exists a strong 
CEF splitting estimated in YbBi2, and we believe that 
increasing pressure can continuously reduce the difference 
in hybridization of the two CEF quasi-quartets. So be- 
low the critical value SV C , YbBi 2 is an insulator with an 
indirect gap as observed in experiments^—, while above 
this critical value SV C this material can transform into 
a heavy electron metal with an enhanced effective mass 
due to the presence of heavy charge carriers. Thus, our 
theory may provide a general microscopic mechanism of 
the pressure induced insulator-to-metal transition in Yb- 
based Kondo insulators/semiconductors. 

Since a constant density of states for the conduction 
electron band was assumed in the above slave-boson 
mean-field calculation, the obtained results are indepen- 
dent of the dimensionality of the model. In order to 
see the special Fermi surface structure of the metallic 
phase, the model Hamiltonian Eq.([T]) is redefined on a 
two-dimensional square lattice system with the conduc- 
tion electron band 

£k = — 2t(cosk x + cos k y ) + At' cosk x cosk y , (13) 

where t denotes the nearest neighbor hopping and t' de- 
notes the next-nearest neighbor hopping. Then the same 
slave-boson mean field calculation can be performed, 
and the insulator-to-metal phase transition also takes 
place for a set of parameters V\ = 0.473, £/ = —0.573, 
t = 0.2 573, and t' = Q.3t when decreasing the parameter 
SV. In the metallic phase, we have calculated the corre- 
sponding Fermi surface structure shown in Fig|3J There 



4 



1.1 
at 

-0.5 
-1.1 



(a) 



-1.0-0.5 0.0 0.5 1.0 

k/ic 



l.o- 

0.5 
0.0 

-0.5 

-1.0- 



(b) 



-1.0-0.5 0.0 0.5 1.0 



1.0 
0.5 
0.0 

-0.5 

-1.0- 



(c) 



-1.0-0.5 0.0 0.5 1.0 



FIG. 3: (color online). Fermi surface structure of the two- 
dimensional model Hamiltonian on a square lattice with Vi = 
0.47), A = 0.173, and ej = -0.57). (a) 5V = 0.087), (b) 
8V = 0.057) and (c) 5V = 0.017). 



exist two Fermi pockets: one electron-like in the center 
of the Brillouin zone and one hole-like in the corners of 
the Brillouin zone. These two Fermi pockets have ex- 
actly the same area in the Brillouin zone. Such a heavy 
electron metal corresponds to a semi-metal. Decreasing 
the hybridization difference 8V below the critical value, 
the sizes of the electron and hole Fermi pockets become 
larger and larger, as displayed in Figj3] 

For the three dimensional Yb-based Kondo insula- 
tors/semiconductors, the Fermi surface of the resulting 
metallic phase should still be given by a hole and an 
electron pockets. As the temperature is lowered enough, 
some instabilities may further appear. Due to the pres- 
ence of the strong mixed valence effect in such systems, 
the additional on-site Coulomb interaction between the 
conduction electrons and localized /-hole should be taken 
into account. In the slave-boson representation, it is 
given by 



(14) 



When the coupling strength Uf c is assumed to be small, 
we can rewrite this additional interaction in terms of the 
hybridized quasiparticles as 



v ki i Mk 2 1 ^k 3 1 Mk 4 1 "k 3 1 "ki l j «k 2 1 



kik 2 k 3 k 4 

+/«k 1 2^k 2 2Mk 3 2t'k42/3k32/ 3 ki2/3k42^k,2 

+ Vki i ^k 2 2 ^k 3 1 Hc 4 2 al 3 ! a kl i Pl 4 2 /3k 2 2 

7 t k 1 27 t k 2 l7 i k 3 27 i k 4 l/5k 3 2^ki2a k4 l Q! k 2 l) , (15) 



where ki + k.2 = IC3 4- k.4 should be satisfied and only the 
two quasiparticle bands, crossing the Fermi energy, have 
been taken into account. a kjl and a k 1 are defined on the 



electron Fermi pocket, while /3 ki 2 and 2 are defined 
on the hole Fermi pocket. Among these residual quasi- 
particle interactions, the first two terms represent the 
intra-pocket scatterings with a small momentum transfer, 
while the last two terms correspond to the inter-pocket 
scatterings with a large momentum transfer. 

According to the recent renormalization group analy- 
sis for a two-band interacting model with electron and 
hole Fermi pockets^, the inter-pocket quasiparticle in- 
teractions will determine the possible instabilities at low 
temperatures. When we set q as a small momentum and 
Q as a large momentum which is the distance between 
the centers of two Fermi pockets, then inter-pocket quasi- 
particle interactions can be approximated as 



-TT (7Wq'lMq2Tiq'2 + V ll P l'X VtfiV^) 
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a q ,i/ 3 Q+q,2^Q +q '2 a qM- 



(16) 



If there is a strong nesting between the hole and elec- 
tron Fermi pockets, this inter-pocket repulsive interac- 
tion will further induce a particle-hole pairing instability, 
corresponding to an orbital-density wave ordering. The 
corresponding order parameter is given by (a; ql /3Q+q2) 
or (/3q +q ' 2 Q q'i)- Such a new type of ordering in heavy 
fermion materials will be discussed in our further inves- 
tigations. 

In conclusion, we have studied the Yb-based Kondo in- 
sulators with a strong CEF splitting in the framework of 
the periodic Anderson lattice model by using the slave- 
boson mean-field approximation. The obtained ground- 
state phase diagram and the indirect gap have demon- 
strated that a second-order insulator-to-metal transition 
occurs via reducing the hybridization difference of the 
two CEF quasi-quartcts. Our theory provides a general 
microscopic mechanism of the pressure induced insulator- 
to-metal transition, because increasing the external pres- 
sure can effectively reduce the anisotropy of the hy- 
bridization strengths of the two CEF quartets experi- 
mentally. The resulting metallic phase has a hole and an 
electron Fermi pockets, which may exhibit an instability 
of an orbital-density wave ordering at low temperatures 
when the inter-pocket quasiparticle residual interactions 
are taken into account. These theoretical results are cer- 
tainly needed to be confirmed experimentally in the fu- 
ture. 
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